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Abstract —In this paper, we study the downlink performance of 
two important 5G network architectures, i.e. massive multiple- 
input multiple-output (M-MIMO) and small-cell denslficatlon. 
We propose a comparative modeling for the two systems, where 
the user and antenna/hase station (BS) locations are distributed 
according to Poisson point processes (PPPs). We then leverage 
both the stochastic geometry results and large-system analytical 
tool to study the SIR distribution and the average Shannon 
and outage rates of each network. By comparing these results, 
we observe that for user-average spectral efficiency, small-cell 
densification is favorable in crowded areas with moderate to high 
user density and massive MIMO with low user density. However, 
small-cell systems outperform M-MIMO in all cases when the 
performance metric is the energy efficiency. The results of this 
paper are useful for the optimal design of practical 5G networks. 

Index Terms —Fifth generation (5G) systems, massive multiple- 
input multiple-output (MIMO), small cell denslficatlon, stochastic 
geometry, large system analysis. 

I. Introduction 

The fifth generation (5G) communications system needs 
to deliver improvement over 4G, e.g. much higher user data 
rate/spectral efficiency, enhanced eco-friendliness and energy 
efficiency, seamless communications, and low latency etc. 11. 
To meet the required performance, densified topologies with 
two main approaches, i.e., massive multiple-input multiple- 
output (M-MIMO) and small-cell densification, are promising 
candidate technologies El. 

In M-MIMO, each base station (BS) employs a large-scale 
antenna array in linear, cylindrical, or other shapes. The large 
number of antennas not only provides more diversity to the 
transmission but also “hardens” the channel, allowing low- 
complexity but sharp beamforming towards the user equip¬ 
ments (UEs) 0. The results are less interference and higher 
spectral and energy efficiencies. The main drawback of M- 
MIMO is the channel estimation/acquisition and feedback, 
which require a high accuracy to attain the gain that M- 
MIMO promises a-ia . Pilot contamination, itself an existing 
problem, also restricts the benefits of M-MIMO ||4|. ITl. 

In contrast, small-cell networks consists of many mi- 
cro/femto cells with limited number of antennas, typically one 
or two, at each access point (AP). Due to the short distance 
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between the APs and UEs, small-cell networks have smaller 
path loss and co-channel interferences while requiring lower 
power consumption, hence improving both spectral and energy 
efficiencies. The gain of small-cell densification, however, 
is subject to more complicated cell planning and inter-cell 
coordination/cooperation 0. 

The evolution of high spectral/energy efficiency for 5G thus 
results in two opposing ideas: concentrating the antennas to 
form M-MIMO BSs, and distributing them to form small cells. 
It is thus natural to ask which approach performs better and 
under which scenarios. Despite many works studying either 
M-MIMO or small-cell networks, there exists few work in 
literature comparing their performance. One of the reasons 
is due to system modeling: no existing framework so far 
allows for an effective collation. In this paper, we aim to 
partially answer the above question using the metrics of 
spectral and energy efficiencies. In particular, we resort to 
stochastic geometry and random matrix theory (RMT) for the 
analysis. A brief introduction of stochastic geometry can also 
be found in 0. 

RMT, also referred as large-system analysis in literature, 
was used to investigate large multi-user linear receivers and 
code division multiple access (CDMA) networks m, m 
and various MIMO setups ina, ini. While exact analyses 
are realizable under finite-dimensional matrix theory m, they 
often lead to complicated expressions hence do not provide 
useful insights. The seminal monograph ifTSll and other works 
have shown that large-system analysis can reveal many more 
interesting properties of communication systems. Furthermore, 
although the analysis is mathematically accurate only for 
large random vectors/matrices, the obtained results remain 
good approximations for finite-dimensional cases. Recently, 
many studies have also advanced the so-called “deterministic 
equivalent” approach to investigate more sophisticated system 

models na. 

M-MIMO, which employs a large number of trans¬ 
mit/receive antennas, is a natural application of large-system 
analysis. A detailed survey can be found in iflTll . which lists 
many key M-MIMO performance analyses considering hard¬ 
ware impairments na, channel state information (CSI) esti¬ 
mation/acquisition 0, 0, and resource allocation/precoding 
schemes ifT^ . 1^ . Note that most of the M-MIMO studies 
only consider a BS-centric analysis in which the relative 
locations of the BSs and UEs are fixed, thus neglecting the 
network topology. In order to compare M-MIMO with small¬ 
cell densification, a proper system modeling is necessary. 

In this paper, we study both small-cell and M-MIMO 
networks using stochastic geometry and large-system tools. 
The aim is to provide an insightful comparison between these 
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two topologies. The contributions of our paper are summarized 
as follows: 

• System Modeling: We propose a model for M-MIMO 
and small-cell networks which incorporates the random¬ 
ness of both transmitting distributed antennas (DAs)/BSs 
and single-antenna UEs. Particularly, the DA/BS and UE 
locations are distributed according to Poisson point pro¬ 
cess (PPPs) with possibly different densities. Therefore, 
there exists a transmission probability that some transmit¬ 
ting antennas/BSs are turned off due to no associated UE. 
We assume a time/frequency division multiplexing access 
(TDMA/EDMA) scheme in which the transmitting nodes 
associated with more than one UE will communicate 
with only one UE in each resource block. Einally, we 
assume non-coordinated DAs/BSs and employ conjugate 
beamforming, which is optimal for non-coordinated M- 
MIMO with single-antenna UEs. 

• M-MIMO and Small-Cell Network Analysis: We study 
the performance of M-MIMO and show that the signal- 
to-interference ratio (SIR) distribution at a typical UE 
is similar to that under a non-fading ad-hoc network. 
Here, the important differences between the two cases 
are the M-MIMO gain and transmission probability. As 
discussed before, the results for fading ad hoc networks 
can be directly extended to the small-cell counterparts, 
with small but important modification again due to the 
transmission probability. We therefore utilize the results 
from l|9] Sections III and IV] to derive the bounds for the 
SIR distribution of M-MIMO and small-cell networks. 
Based on these analyses, we further obtain the bounds for 
both the average Shannon and outage rates of the two sys¬ 
tems. To better characterize the small interference regime, 
the performance of M-MIMO and small-cell networks is 
also derived under the assumption of asymptotically small 
UE density. 

• Comparison between M-MIMO and Small-Cell Den- 
siflcation: We reveal that small-cell densification yields 
better rates than M-MIMO when the user density is 
moderate or large compared to the transmitting node 
density and number of antennas. However, M-MIMO 
outperforms small-cell systems under the asymptotically 
small UE density regime. This implies that there exists a 
UE density threshold lower than which we should employ 
M-MIMO and higher than which small-cell densihcation 
is more preferable, when spectral efficiency is the perfor¬ 
mance metric. We further compare the energy efficiency 
performance and show that small-cell always outperforms 
M-MIMO. The results are therefore useful for the optimal 
design of the upcoming 5G heterogeneous networks. 

The rest of the paper is organized as follows. Section HI] 
describes our M-MIMO and small-cell network models, which 
incorporate the randomness of both transmitting nodes and 
UEs, for a fair comparison between M-MIMO and small-cell 
densification. In Sections Hill and HY] we investigate the M- 
MIMO and small-cell systems introduced in Section HD based 
on the results from ^ Sections III and IV]. The spectral and 
energy efficiency comparisons between M-MIMO and small¬ 


cell densification are given in Sections IVl and IVll respectively. 
Einally, Section rvHI concludes the paper with some remarks. 

Notations: Scalars and vectors/matrices are denoted by 
lower-case and bold-face lower-case/upper-case letters, re¬ 
spectively. The conjugate, transpose, and conjugate transpose 
operators are denoted by (•)*, (-j^, and {■)^, respectively. 
\A\ij stands for the {i,j)th element of the matrix A. E^]-] 
denotes the statistical expectation over a random variable X. 
Tr(-)’ II ■ lli^’ det(-) represent the trace, Erobenius norm, 
and determinant of a matrix, respectively. Im denotes an M- 
hy-M identity matrix. and denote the space of 

x-hy-y real and complex matrices, respectively. Einally, the 
circularly symmetric complex Gaussian distribution with mean 
/i and variance is represented by CJ\f 

II. Stochastic Modeling for Massive MIMO and 
Small-Cell Systems 

A. M-MIMO System 

We consider a M-MIMO system where BS and UE locations 
are distributed according to homogeneous PPPs $c',d and 
of intensities Ac,b = Af, and Ac,ii = ^d,u = A„, respectively. 
Each UE is associated with the nearest BS each with M 
antennas. Here, we assume that each BS is scheduled to serve 
only one UE at each time/frequency slot due to time/frequency 
division multiple access (TDMA/EDMA). In each slot, the 
received signal of the UE associated with BS b is given as 

yC,b = hQijl^Xc,b + + n-C,bj (1) 

6'GB/{h} 

where B is the set of transmitting BSs; ^ is the channel 
between BS b' and the associated UE of BS 6; and xc.fe' is the 
transmit signal from BS b'. Here, hc,b',b = \/^c,b',bhc,b',b'^ 
^C,b',b = d^^b' b i^ n large scale fading including path loss 
with ji denoting the path loss exponent; hc,b',b follows a 
statistical distribution hc,b',b ^ CA/'(0,1), \/b,b' G B\ and 
nc,b is the additive white Gaussian noise (AWGN) distributed 
as CA/'(0,cr^). Eor a fair comparison with small-cell systems, 
we assume that each BS has a total power constraint of MPt 
and therefore = MPt- 

We consider a practical regime in which each BS only 
knows the channel between itself and the associated UE via, 
e.g., channel estimation and feedback. Under such assumption, 
each BS implements the conjugate beamforming to transmit 
signals to its respective UEs. 

Eig. m shows one channel realization of the M-MIMO 
system model. Here, the service area is assumed to be a square 
region of 20 x 20. The BS and user densities are Ac,b = 0.05 
and A„ = 0.15, respectively. Note that as the UE density is 
not sufficiently larger than the transmitter density, some BSs 
are not associated to any user. The corresponding cells are 
marked by green color for a clear illustration. 

B. Small-Cell System 

Eor a fair comparison, we consider a small-cell system 
where distributed antenna (DA) and UE locations are dis¬ 
tributed according to homogeneous PPPs ^D,b and „ of 
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Fig. 1: A channel realization for either M-MIMO or small-cell 
systems. Here “o”, “x”, and straight lines denote the transmitter 
(BS/DA), user, and boundaries between cells. The cells without any 
user are marked by green color. 

intensities X^.b = MXb and ^ = A„, respectively, in the 
Euclidean plane. Each DA is associated with the nearest UE. 
Eor simple description, we consider one resource block (e.g., 
frequency and time), which is taken by only one UE, i.e., one 
DA serves one user at a time. The received signal at user k is 

yD,k = hD,k,kXD,k + ^ hD,k,k'XD,k' + Uu^k, ( 2 ) 

fe'G/C/{fe} 

where 1C is the user set, hu k,k' is channel gain between 
user k and selected DA for user fc'; XD,k is the transmit 
signal from user fc; and nu^k is the AWGN CA/'(0,cr^). 

Here, hD,k,k' = \/Cn,fe,fe'^n,fc,fc'; ^D,k,k' = djj^k,k' ^ 

scale fading including path loss with /r denoting the path 
loss exponent; and hD k,k’ follows a statistical distribution 
hD.k,k' ^ CA/'(0,1), Vfc, fc' G JC,k' ^ k. We assume equal 
power constraint at each DA, i.e., W\XD,k\'^ = Pt- Note 
that we can consider Eig. [T] as a channel realization for the 
small-cell system model with the DA and UE densities being 
Xdx = 0.05 and A„ = 0.15. 


111. Massive MIMO System Analysis 

In this section, we investigate the performance of M-MIMO 
systems based on the results in our paper ||9l Section II]. Eirst, 
we derive the distribution of the signal-to-interference-plus- 
noise ratio (SINR) and SIR. Note that if the UE intensity A„ is 
not sufficiently larger than Af,, then some M-MIMO BSs may 
not possibly be associated to any UE, and thus, do not transmit 
signals. Therefore, it is necessary to study the transmission 
probability of each BS. The following result provides the area 
distribution of Voronoi cells. 

Lemma 1 ( KITil): Given a PPP $ with homogeneous den¬ 
sity A. The probability distribution function (PDE) of S, 
the area of a typical Voronoi cell formed from <1), can be 
approximated as 

r(3.5) 


By using Lemma [T] the transmission probability of BSs is 
approximated as ll^ 

Note that given the transmission probability ec, the trans¬ 
mitting BS locations are thus distributed according to a PPP 
$c,f) with density Xc,b = f-cXb- 


A. Conjugate Beamforming - Large System Analysis 

Since all BSs employ conjugate beamforming, the transmit 
signal xc,b is 


\/^T^hc,b,b 

Xc,b = Wc,bSC,b = -fn:-f]- sc,b, 


( 5 ) 


\\hc,b,b 

where wc,b — \/MPrhc ,b,b /\\hc ,b,b\\ C is the con¬ 

jugate beamformer of BS b. As a consequence, the received 
signal of the UE associated to BS b is given as follows 


yC,b = sj MPT\\hc,b,b\\^C,b 

^ , ,,, 

+ -nr-n-+ nc,b, (6) 

The SINR of the UE associated to BS b therefore can be 
expressed as 

MPT\\hc,b,bV 


SINRb = 


itrn ^C,b',b^c,b',b'hc^b',b'^C,b',b 2 

2^ M Ft - T7-. - T77. -h cr„ 


b'eB/{b} 




^C,b,b 


he, 


b,b 






E 

b'eB/{b} 


j-ti ^C,b'.6^C,6'.b'^C.h'.b'^C.b',b , 
^C,b',b T + 


hc,b',b' 


MPt 

( 7 ) 


Note that the vector series 


he 


s/M 


e C 


Mxl 


satisfy the 


condition of iflhl Theorem 3.4]. Therefore, we have 

2 


hc,b,b 




( 8 ) 


Based on ([8]), the SINR of the UE associated to BS b is 
expressed as 

Qm = SINRb 

2 


^C,b,b 


hc,b,b 


/M 


E 

b'eB/{b} 


^C,b',b 


hc,b',bhc,b' ,b' h^ k',b'hc,b' ,b 


M—foo 


M 

Mdclb 


hc,b',b' 


M^Pt 


,•2 ’ 


( 9 ) 


fe'eSc.d/{b} 


where hb' = ; hbds are i.i.d. random variables 


each distributing as CAf (0^ 1). The form @ is similar to that 
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investigated in |[9l Section V]. It has been shown that 
is difficult to study, in particular to derive the cumulative 
distribution probability (CDF) bounds. As discussed in El 
Section V], we can well approximate 




C,b',b- 


b'&<i>c.d/{b} 

As a consequence, we have 


b'e$C,d/{6} 


Qm 




E 

6'6$c,d/{6} 


^C,b',b 


MPt 


( 10 ) 


( 11 ) 


B. Probability Distribution 

We note that the form of the SINR in (fTTT l is similar to 
that investigated in ® Section II]. Therefore, we can apply 
the results for the outage probability, average achievable rate, 
and outage rate in ® Section II] here. The only but important 
difference is that the interfering BSs’ locations are distributed 
according to the PPP with density Ad, while the PDF of 
the nearest distance follows that from PPP with density 
Afc. The following corollaries follow same reasonings as in ||9l 
Theorem 1, (9), and (10)]. The proofs are given in Appendix 
lAl where we only highlight the differences for brevity. 

Corollary 1: Assume that the transmitting BSs follows a 
PPP <hc,d with density \c,d- The Laplace transform of 
can be expressed as 


C 


Q-i (s) = TrAfa 


exp 


M^Pt 






Corollary 2: Given that M —^ oo, the lower and upper 
bounds of the SIR distribution are given as follows 


Fhlid) 


0, 0<(7<M,(13a) 

1-^ ^ (13b) 

ecM r-qr- + 1 — ec 




1 

1 + ’ 


(14) 


Similar to El (10)], (O is rigorously proved only for: (a) 
q > M, and (h) q < M with g —>■ 0 or g —M. The following 
result is a counterpart of ||9l Corollary 1] for massive MIMO 
when M ^ oo. The proof can be deduced from Appendix 
lA-BI and thus is omitted for brevity. Note that Corollary [3 
does not require the approximation (fTOl i since [MSIR]]"^] 
already averages over the channel fading. 

Corollary 3: Assuming a very large M, the expectation of 
MSIR^" is given as E<j. [MSIR,^^] = 

Remark 1: It is important to remind that Corollaries|2]and[3 
are derived by assuming that the noise power is negligible, i.e., 
— 0, and considering the SIR instead of SINR. In small 
UE density network, i.e., A{, ^ A„, such assumption is not 
accurate anymore since the interfering sources are far away 
and the interfering power is small. The noise effect therefore 
is more important. As a consequence. Corollaries |2] and [2 are 



Fig. 2: Simulated SINR CDF and its bounds ([T3), ([T4j for M-MIMO 
with M = 64, p = 3.7, Af, = 1, and Pt/ct^ = 15 dB. 



Fig. 3: Average of MSINR~^ for M-MIMO with M = 64, p = 3.7, 
A6 = 1, and Pt/ct^ = 15 dB. 


not accurate in UE-sparse networks. Eurther results for M- 
MIMO with asymptotically small UE density will be presented 
in Section UlI-DI 

In Eigs. |2] and 12 we confirm the validity of Corollaries 
12 and 12 respectively. Erom Eig. [2 it is observed that the 
assumption that the transmitting BSs follow a PPP ^c,d with 
density Xc,d = ec^b with ec given in (|4li is quite accurate. 
However, the upper bound is only close to the exact SINR 
distribution when A„ is sufficiently large {Xu > Xb by intensive 
simulations), as shown in Pig. |2 As Xu decreases compared 
to Xb, the upper bound (fT4l i becomes looser. The reason is that 
some inequalities, used to derive the upper bound in Appendix 
EH] are not strict especially when ec is small. The lower 
bound, however, is not affected by the UE density due to its 
simple derivation. 

C. Average Achievable and Outage Rate 

Based on Corollary |2 we investigate the average Shannon 
and outage rates for M-MIMO systems. We note that due 
to TDMA/PDMA, each user is only served for a fraction of 
time or a sub-band. Assuming that all users have the same 
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transmission priority, that fraction can be defined as 
where is the number of users in cell b. The following 
result establishes E 

Lemma 2: The expectation of the number of users in a 
typical cell h is given as 

E[Nt,u] = ^ = ^. (15) 

Proof: Please refer to Appendix IbI ■ 

As expected, the average number of users in a cell is 
proportional to the UE density while inversely proportional 
to the BS density. The simple linear scaling is surprising, but 
intuitive. The reason is due to the homogeneity of both BS and 
UE PPPs, which does not create clusters of either BSs or UEs. 
Therefore, with A„ UEs and Af, BSs in a region of 1 km^, we 
expect to see the UEs homogeneously distributed across all 
BSs, which leads to A„/Ab. By comparing the simulated and 
analytical values of E [A"f, „], intensive simulations confirm the 
validity of Lemma |2] The comparison is not presented due to 
the space constraint. 

Note that the rate of an UE associated to BS h is given as 

Rb,u = log 2 (1 + SlNRf,), (16) 

Plb,u 

where the fraction -rf— is due to TDMA/EDMA. We thus 

^b,u 

obtain the average Shannon rate of a typical user as follows 


(a) 


Rm = = IE 


N, 


b,u 


E [log 2 (1 + SlNRb)] . 

(17) 


Here, (a) is due to TDMA/EDMA and the fact that the 
random variables SINR;, and A"f, „ are independent given that 
TDMA/EDMA is applied. 

Based on Corollary |2] and Lemma |2] we can derive the 
bounds for the average achievable rate and the outage rate as 
follows. The proofs are similar to that of ||9] Lemmas 3 and 
4] and omitted for brevity. 

Corollary 4: The average user rate Rm is bounded above 
and below by R}^ and 7?^®, respectively, where 

i?™=min|l,^| ^log2(l + M) 

/og2(i+M) — 1)2/^ + 1 — ec 



„iR . dt 

l + ec/3n-fdM-2/M(2*-l)2/M' 

(19) 

Corollary 5: The outage user rate ORMiv) is upper- and 
lower-bounded by OR^® ( 77 ) and ORj)f ( 77 ), respectively, where 


min 



l0g2(l + 7 ), 

l0g2(l - 


v) 


ec{ 


MJ 


1 - ec 


0 < t] < M, (20a) 
M < 77 . ( 20 b) 



Fig. 4: Average user rate and its bounds GD and (Ell for M-MIMO 
with respect to M. Here, p = 3.7, Xb = 1, and Pt jrP' = 15 dB. 


0 R“ ( 77 ) = min 


^ ^ 1 log2 (1 + ??) 

l + ec^n-fd(l)'/'^' 


( 21 ) 


Here, the term min|l,-^| in Corollaries |4] and |3 is due 
to the fact that TDMA is applied only when A^f, „ > 1, and 
we have approximated 


E 


1 

d^b,u 


1 _^b 

IE [TVfe.u] A„ 


( 22 ) 


In Eigs. |4] and |5] we compare the average and outage user 
rates with the bounds given in ( fT8]) -(l2Tll. We first observe that 
the scaling laws of the bounds are closely matched with the 
true rates. The gap, however, seems to increase as the UE-BS 
density ratio becomes smaller. Nevertheless, the bounds are 
useful for further studies on the M-MIMO performance under 
various network models. 


D. Asymptotically Small UE Density 

In this section, we investigate the performance of M-MIMO 
under the asymptotically small UE density regime. As the 
transmission probability ec approaches 0, the interference also 
goes to 0 since the interfering BSs are farther apart from each 
other. In such cases, we can ignore the interference but not 
the noise due to its dominant effect. The SINR can now be 
approximated as the SNR and expressed as 


Qm = SINRb « SNRfa = 


MUrdc.&.f) ^c,b,b 


a.s.^ 


(23) 


Since the PDE of dc,b,b is 

fDc,b,b(x) = 27rAba7e“’"^'’^^ (24) 
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Fig. 5: Outage rate and its bounds l l20t and l l21t for M-MIMO with Fig. 6: Simulated SINR distribution, analytical result Oil ), and 
M = 64, jj, — 3.7, Xb = 1, and Pt/i^^ = 15 dB. bounds j33b . l l34b for small-cell system with Xofi = MXh = 64, 

= 3.7, and Pt/ct^ = 15 dB. 


it is straightforward to derive the SINR distribution as follows 
Pr{QM < 9 } = Pr < dc,b,b > 


\ 9Crn / J 


= exp —irXb 




(25) 


\ J 

Therefore, the average Shannon and outage rates are given 


as 


i?M — 


1 — exp -nXb- 


Cr4/M(2t _ 1)2//.^ 


dt, (26) 


ORM(f?) = log2(l + v) 


1 — exp —irXb 


^ 4 /Mp 2 /M\ 


(27) 


IV. Small-Cell System Analysis 


Similar to M-MIMO, if the UE intensity Xu is not suffi¬ 
ciently larger than Xb, then some DAs in the small-cell system 
may not be associated to any UE, and thus, do not transmit 
signals. Applying Lemma [T] we can derive the transmission 
probability of a DA in our system as 1221 




1 - 



A. 

3.5MXb J 


(28) 


The transmitting DAs thus can be modeled as a homoge¬ 
neous PPP b obtained by thinning the DA PPP ^D,b with 
a probability eo that DA is activated. In this case, the SINR 
of a typical UE k can be described as 


Qsm — 


\hD 


k.k 


\^dZPT 




k'£^D,b/{k} 


lk,k' 


PTdZ'+< 


(29) 


Using CD, we can obtain the coverage probability of a 
typical user k in Lemma l3] ll2Tl . 

Lemma 3 ( Given that the transmitting DAs are 

modeled as PPP & and both desired and interfering signals 


obey Rayleigh distribution, the coverage probability for user 
k is expressed as 


Pr (Qsm > q) 


irMXb 

irMXbV 


exp <; - 

Pt 



(30) 


A direct corollary of Lemma 6 is the SIR distribution for 
the case a'^jPr = 0 given as follows. 

Corollary 6: The SIR distribution of a small cell system is 
given as 


PQsMiq) = 1 - 


1 

1 + eoq^/^p 


(31) 


Again, we obtain the SIR CDE, average achievable rate, and 
outage rate bounds for the small-cell system. The differences 
between the following Corollaries [T] [8] and Lemma 6] are 
that the interfering cells follows a PPP with density euAb, and 
each DA might need to serve multiple UEs via TDMA/EDMA. 
We note that the expectation of the number of UEs associated 
with a typical DA k is given as (cf. Lemma [JJ 




^D,u 

)^D,b 


MXb 


(32) 


The proofs for the following results are similar to those of 
||9l Lemmas 5 and 6], and thus omitted for brevity. 

Corollary 7: The lower and upper bounds for the distribu¬ 
tion of Qsm defined in (l23 l are given as 




;iUB 


( 9 ) = 1 - 


, , 2 £Dg^p (2 

1 


where \h\’^ is exponentially distributed. 



(33) 

(34) 
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Fig. 7: Average user rate and its lower bound l l35t with respect to Fig. 8: Outage rate and its lower bound ([36} for small-cell system 

M for small-cell system. Here, \h = 1, = 3.7, PtIo'^ = 15 dB, with A_d,6 = M\h = 16, jJ, — 3.7, and PtI^^ = 15 dB. 

and \o,h = MA(,. 


Corollary 8: The average Shannon rate i?sM(ct) and the 
outage rate ORsmCck, v) are bounded below, respectively, by 


-^SM = min 

ORsm(' 7) = min 


MAfc ] ^ 

) Jo l + e_D/lfd(2* 

MXb \ log2 (1 + t?) 

Au J 1 + 


(35) 

(36) 


Again, the term min 1 1 ,| is due to the fact that 
TDMA/FDMA is applied only when Nk^u > 1, and we have 
approximated 


1 _ 1 _ MXb 

Nk,u E [Nk,u] Xu 


(37) 


In Fig. |6l we compare the SINR distribution, analytical 
result (OTI) . and bounds ([33}, (134} for two small-cell systems. 
It is observed that the upper-bound (l3^ . although simple, is 
quite tight when compared with the lower-bound ([33}. Figs. 
[7| and [S] illustrate the average and outage user rate for various 
small-cell systems. We also show the analytical rate lower 
bounds given in Corollary [8] for comparison. In contrast to the 
M-MIMO case, the bounds for small-cell systems are quite 
tight. Note that we do not derive and present the analytical 
upper bounds for the rates here since they are complicated 
and of limited use for subsequent analyses. 


A. Asymptotically Small UE Density 

Under the asymptotically small UE density regime, the 
interference is negligible and the SINR of small-cell system 
can be expressed as 

\hD,k,k\ dkk^T 

Qsm = SINRfe Ri SNRfe = J—. (38) 

<^n 

Since the PDF of dk,k is 

fD,,, (x) = 2TrMXbxe -^^^'>^", (39) 


it is straightforward to derive the SINR distribution as follows 


Pr{(5sM < q} = E/j 


exp -TrMXb- 




(40) 


\ (Tn q 

where the random variable h ^ CAf{0, 1). 

Therefore, the average Shannon and outage rate are 


Rsm — 


1-E?, 


32/m 


exp --KMXb—^j 






dt, 
(41) 


ORsm(?7) 

= l 0 g 2 (l + 77 ) ( 1 - Eft, 


32 /M 


exp I —irMXb — 


\h\y^ 




(42) 

To derive the bounds for (|40}- (|42} . we note that the function 
exp(—aa;^) with /3 < 1 is convex in x > 0. Therefore, 
applying Jensen’s inequality, we obtain 

^ 2 //i-n 


> exp --KMXb ^ ^ ^ 

v Ori q^^^ 


pj/Mr(l + g)' 

n^l M 


= exp I —TiMXb 


= ^qL(7)- 


(43) 

The upper bounds for i?sM and ORsm( 7?) are therefore given 

p2/Mr(i +1) 


pUB _ 
TtSM — 


1 — exp -irMXb- 


f^4/M(2t_l)2/M^ 


dt, 

(44) 


OR^m(7?) 

/ / 4- 

= log 2 (l + 77 ) I 1 - exp ( -nMXb- ^ 




(45) 
























































V. Rate Comparison for M-MIMO and Small-Cell 
Systems 

M-MIMO and small-cell represent two densification ap¬ 
proaches in 5G communication systems. The lack of a perfor¬ 
mance comparison between them is even more surprising given 
the fact that each approach has attracted significant attention 
in the literature. In this section, we exploit the newly obtained 
results in Sections m and HV] to draw several observations on 
M-MIMO and small-cell networks, revealing which approach 
is better and under which setup. Here, the metric is the outage 
rate for simple arguments, but our observations can be readily 
extended to the average Shannon rate. The key parameter here 
is the UE density A„. 


A. Very Large UE Density: > MXh 

In this case, the lower bounds of the outage rate are reduced 
to 




Xb log 2 (1 + t?) 

1 -I- ec/3n-fdM“2/M,y2//i ’ 


(46) 


ORZiv) = 


MXb log 2 (1 + v) 

1 -I- ’ 


(47) 


where ec ~ 1 ~ ~ 1- Simple manipulation 

shows that ORj^ iv) > ORiif (.v) iff 

M + > 1 + (48) 


As a loose approximation, (l48l l holds if > 

eD/?fd. which is true since /3n-fd > Pfd, ec ~ 1 > and M is 
large. We thus observe that a small-cell system outperforms a 
M-MIMO counterpart when the UE density is very large. This 
is due to the effect of multiplexing large number of users under 
TDMA/EDMA, which diminishes the interference-mitigation 
benefit of M-MIMO. 


B. Intermediate UE density: Xb < Xu < MXb 


Since 


A„ 


3.5MA|, 

CD ~ 1 - ( 1 


is small, we have 

Ai, 


-3.5 


Aii 


MXb' 


(49) 


3.5MXb, 

Based on (l2Tl l and (|36] |. we can express the outage rates as 
Xb log2 (1 + ry) 


ORm (^) = 


Au 1 + ecPn-fdM ’ 


(50) 


ORS(^) = 


l0g2 (1 + ??) 


1 -f en/ifdp^/^ 
The condition for ORsm(? 7 ) > ORj^f (p) is 


(51) 


(52) 


We note that (|52] | is true since 1 > -^ and M > 
Therefore, s small-cell system outperforms a M-MIMO coun¬ 
terparts when the UE density is moderate. Similar to Section 



Fig. 9: Comparison of the achievable user rate between M-MIMO 
and small-cell systems with A^ = 1. 


IV-AI the reason is again due to the effect of multiplexing a 
large number of users under TDMA/EDMA M-MIMO. 

Remark 2: In Sections [V- Al and IV-Bl a stronger and stricter 
result which states that ORs^(77) > OR™(p) is also provable 
using the same arguments. However, we have used the lower 
bounds ORjjf (ry) instead due to the symmetry and for the ease 
of explanation. 


C. Asymptotically Small UE Density: Xu A;, 

Assuming a small UE density, it is difficult to compare 
M-MIMO and small-cell systems since the M-MIMO lower 
bounds are not close to the exact distribution or rates. In this 
subsection, we consider the asymptotic regime where Xu is 
very small compared to Xb to reveal some insights for the 
small UE density case. Note that 


exp 


T:XbPy^ 



< exp 




MT 



(53) 


when p, < 4 and M is large enough, i.e., M^ > T ^1 -F 
Erom ( l26b . ( l27b . (l44l) . and (1451 ). it is straightforward to prove 
that Rm > RZ and ORm(p) > ORZiv)- 

The above result essentially asserts that M-MIMO outper¬ 
forms small-cell densification, albeit when the UE density is 
asymptotically small and p < 4. Combined with observations 
from Sections IV-AI and IV-Bl we note that there exists a 
threshold with smaller UE density than which we should 
employ M-MIMO and larger UE density than which small-cell 
densification is more preferable, in terms of spectral efficiency. 
As an example, we compare the achievable user rates of M- 
MIMO and small-cell systems in Eig. |9] Note that as the path- 
loss exponent p approaches 4, Xu needs to be smaller so that 
M-MIMO rate can outperform small-cell counterpart. When 
p > 4, the M-MIMO rate is always worse than that of small¬ 
cell, as suggested by (l5Tt . 

We would like to highlight here that our study have not 
considered the deployment cost, the signaling overhead, the 
handover implementation, etc. We have also not considered 
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multiuser beamforming in M-MIMO. Due to the space con¬ 
straint, we leave such interesting and important investigations 
as our future works. 

Note that we can consider small-cell and M-MIMO systems 
as two extrema of a balancing problem, where the number 
of BS antennas and BS density can vary while their product 
always equals to MXb- In other words, the problem is to 
distribute the antennas from a pool with density MA& so 
that the resulting rate performance is optimal. The analysis 
of this problem is partially covered by Sections [HI] and |IV] 
The above arguments suggest that the optimal setup should 
be a moderate M-MIMO system with sufficient beamforming 
gain and a sufficiently dense BS deployment. 


VI. Energy Efficiency Comparison for M-MIMO 
AND Small-Cell Systems 


Apart from spectral efficiency, energy efficiency (EE) is also 
an important factor for evaluating future communications net¬ 
works due to several environmental concerns. In this section, 
we compare M-MIMO and small-cell systems using EE as the 
metric. In the literature, a conventional definition for EE is 


d- A ^AP/BS 


(54) 


where Pap/bs Pap/bs are the transmit power and sum- 
rate of the corresponding AP/BS to all users, respectively. 
However, for ease of analysis, we consider the following EE 
definition based on the outage rate 


EE(ry) 4 

Pap/bs 


(55) 


We can interpret EE(? 7 ) as the energy efficiency of the 
AP/BS given that the user SINK satisfies the constraint SINK 
> f]. Observations for EE based on dSli can be similarly 
obtained. Note that due to the transmission probability, the 
expected transmission power for each M-MIMO BS and small¬ 
cell AP are given as ecMPr and edPt, respectively. We 
first consider the moderate to large UE density regime where 
Xb < Xu and obtain 


EE^m(^) = 


ORsS(^) 

cdPt 


OR^ OR^ 

^dPt ecMPx 


EE“(77), 

(56) 


since eo < ec, M ^ 1, and OR^^{r]) > ORj;f(p) from 
Sections IV-AI and IV-BI 

Eurthermore, under the asymptotically small UE density 
regime with A„ <C Af,, we have ec ~ and ~ 'M\~- 
Therefore, ec ^ sd- Erom ( |26] |, (l27l i, djlT i. and (l42T i. it is 
straightforward to show that EEsm > EEm and EEsm(i 7) > 
EEM(ry) by noting that the exponential function approaches 0 
as M grows large. 

The above result reveals that small-cell systems are more 
energy-efficient than M-MlMO counterparts with large, mod¬ 
erate, or asymptotically small UE densities. Strictly speaking, 
such results do not cover the small UE density regime. 
However, we conjecture that small cells outperforms M-MIMO 
in terms of EE under all cases. This conclusion is likely to 


hold even with multiuser beamforming/precoding applied at 
each M-MIMO BS. 

Erom Sections |V] and |Vll we observe another important 
aspect of the trade-off between M-MIMO and small-cell den- 
sification. Particularly, even through M-MIMO have a higher 
spectral efficiency as under small UE density regime, small¬ 
cell densification might still be a preferable option for certain 
communication systems due to its higher EE performance. 

VH. Conclusions 

In this paper, we have compared the spectral and energy 
efficiencies of massive MIMO and small-cell systems. Par¬ 
ticularly, we have derived SIR distribution bounds for both 
systems, based on which the average Shannon and outage rate 
bounds are obtained. We have also analyzed the performance 
of M-MIMO and small-cell under asymptotically small UE 
density regime, which represents UE-sparse networks. The M- 
MIMO and small-cell systems were then compared in terms 
of spectral and energy efficiencies. Eor the rate performance, 
we observe that M-MIMO surpasses small-cell densification 
when the UE density is small compared with BS/AP density 
and number of antennas, and vice versa. However, small¬ 
cell network yields better energy efficiency than M-MIMO 
counterpart under all cases. The results of this paper are useful 
for the optimal design of practical 5G networks and other 5G 
system performance analyses. 


Appendix A 

Proof of Corollaries [T|and|2] 

A. Proof of Corollary |7] 

Given that dc,b,b = x. The Laplace transform is given as 

E (26)] 

Cq~i {s\dc,b,b = x) = exp I - 7rAc.&a;^(l - 



(57) 


where 6 = Now we observe that dc,b,b is the distance 

from an arbitrary origin to the nearest point of a PPP with 
density Xb- The PDE of dc,b,b is expressed as follows 


= 27rAba;e 

The Laplace transform of SINR]]"^ 


is thus obtained as 


(58) 


(s) 


pOO 

/ Esinr-1 is\dc.b,b = x) fDc,b,b ix)dx 
Jo 


= 27rAb / xexp<{ - —x>^ 5 - t:Xc, bX^ 


1-b 


Xb 

^C,b 


- 




= 7rA& 


i-- 


M^Pt 


1 + 


ec 


— e 


-s/M 




This concludes the proof of Corollary [T] 
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B. Proof of Corollary \2} 

For ease of presentation, we will consider Qm/M. The 
lower bound is given as 

Pgl/AiiQ) = / = x)fDc,,,d^)dx 

Jo 

(l - 2TTXbXe-^^^^"dx 

^ (60) 


10 
= 1 - 


g^g2//x _|_ 1 


FQ^/Miq\dc.b,i, = x)<l- 


p — 2 

We therefore obtain 

Fq^/Miq) 

< 2TrXb (l - xe-^^^^^dx 


(2TrXb)^ece^>^-^q^/>^ 


H — 2 


= 1 - 


1 


ec(Cg)^/'" + 1 - ec 


1 


(jO _ 


^ ^ (ec(C<?)^/^ + 1 - ec)" 

1 


(fcq^l'^ + 1 - ec) 

2^2/M-lgg,g2/;. 1 


(/.-2)e-2 


B 2 ^4/^1 _|_ 2 ^ _ 

1 


(/r - 2)fi+2/M ecg^/^ + 1 - ec ’ 
where in (a) we have defined 

+ 1 - ec. 


C = 


ecq^/^ + 1 — ec 


fit) 


(P - 2)g - 2 
(/r - 2)fi+2//^ 

f is ^ 

guidance for our subsequent analysis. 


Specifically, from (l62l i. we will show that 
1 ^ 1 


1 - 


> 1 - 


1 + ec(Cog)^/'^ + 1 - ec 


d 2 (ec(Co<7)^/^ + 1 - ec)" 


where = ^0^- After some manipulations, it is observed 
that ( |66] ) holds true if and only if 


Now assume that dc,b,b = x, q > 1, and ^ > 1. Similar 
to (39)], the distribution for massive MIMO case can be 
bounded as 


2 ^ 

^Jj Zd ^ z 

< (/5n-fd - q^^*" + 1 - ec + icifoqf'^^ 

+ (l-ec)(/3n-fd-?o^'‘)g'/^ (67) 

Now note that with /3„-fd = = 0^, we 

have 

^n-fd - = —^/3n-fd, (68) 

fl Z qo jl -r z 


which leads to 


, (61) 




Furthermore, we obtain 


^g'^/^ec/Jn-fd- 


p- 2^° 


(69) 


ec{^oq0>^ + (1 - ec) (/3„.fd - fo^'^) q^'^ - 


~ '^'^''^^/r^- 2 ) 2 /M + 2) + (1 - ec) 


2(/r + 2) 


- 2 


(70) 


(62) 


(63) 


(64) 


and (&) is due to the fact that ^2 /m-i < ^ 2 //i-i ^ > 2 

and 0 < ec < 1. Note that ^ has the same support as 
i.e., ^ G [l,c»). The next step is to find the maximum of the 
function 


(«c>. + (l-cc)^)>0. 

Combining (l62 l and (iTOl) we thus showed that (l67l) holds 
true. Therefore, for g > 1, we conclude that 

Fq^/miq) < 1 - 1 + ecPn-uq^^^^ ^ ^Qm/m(7)- (71) 

Similar to 12 Lemma 2], we can show that 

Fq^/Miq) < 1 - ——- -f - jj- (72) 

1 + ecPn-fdr/'" 

as g —!► 0 or g < 1 but g —1. Intensive simulations show that 
(I72I 1 holds true for all intermediate values 0 < g < 1 as well. 
This concludes the proof of Corollary |2] 

Appendix B 
Proof of Lemma|2] 

Applying Lemma [T] the PDF of S, the area of a typical 
Voronoi cell formed from $c,&, can be given as m (12)1 


o c:3.5 

\3.5 2.5 -3.5A6X 

r(3.5)^'’ ^ " 


(73) 


(65) 


Since the UEs are distributed according to a homogeneous 
PPP with density Xu, the probability of the number of UEs in 
an area x is given by 


given that ^ > 1. The optimal is This value is a 


Pr{A^„ = njarea = a;} = e 


\n^r, 
— XuX '^u'^ 


nl 


I ■ 


(74) 


( 66 ) 


The distribution of the number of UEs is thus equal to 

pOO 

Pr{A^b,u = ^} = / ^r{Nb,u = n|area ^ x}fs{x)dx 
Jo 

q ir3.5 \n poo 

= (75) 

r(3.5) n! Jq 
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Note that by using Il24l (3.371)], we can express the distribu¬ 
tion of the number of UEs as 


Pr{iVb,u = n} 

^ (2n + 5)!! 3.53-5 „ / 1 

n! 15 X 2’" ^ -f 3.5Ab 


n+3.5 

(76) 


This expression, however, is difficult to manipulate. In the 
following derivation, we use (fTSl l instead. The expectation of 
is given as 


E [A^b.u] = ^ nPr{Nb^u = n} 

n—1 


-A3-5 


r(3.5) ,,0 

W 3.5^'® 3 5 


\n—l 
poc 


n\ 


^2.5^-iXu+3.5\,)x^^ 


r(3 


-A^^A^I (77) 
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where (a) comes from the fact that (^) 

is due to ll24l (3.371)]. This concludes the proof of Lemma|2] 
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